A Survey of 286 Virgo Cluster Galaxies at Optical griz and 
Near-IR iif-band: Surface Brightness Profiles and Bulge-Disk 

Decompositions 



Michael McDonald 

^ Department of Astronomy, University of Maryland, College Park, MD 

O 

CN mcdonaldOastro . umd . edu 

, ^ Stephane Courteau 

^ Department of Physics, Engineering Physics and Astronomy, Queen's University, Kingston, 

ON, Canada 

^ courteauOastro . queensu . ca 

^ R. Brent TuUy 

O Institute for Astronomy, University of Hawaii, 2680 Woodlawn Drive, Honolulu, HI 

c3 tully@ifa.hawaii.edu 



Joel Roediger 



C\| and 
> 
<N 
O 
O 

O Department of Physics, Engineering Physics and Astronomy, Queen's University, Kingston, 

ON, Canada 

O 

^ j roedigerQastro. queensu. ca 

jxl ABSTRACT 

We present g,r,i,z and H-band surface brightness profiles and bulge-disk de- 
compositions for a morphologically-broad sample of 286 Virgo cluster catalog 
(VCC) galaxies. The H-band data come from a variety of sources including our 
survey of 171 VCC galaxies at the UH 2.2-m, CFHT and UKIRT telescopes, 
and another 115 galaxies from the Two-Micron All-Sky Survey (2MASS) and 
GOLDMine archives. The optical data for all 286 VCC galaxies were extracted 
from Sloan Digital Sky Survey (SDSS) images. The H-band and SDSS griz 
data were analyzed in a homogeneous manner using our own software, yielding 
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a consistent set of deep, multi-band surface brightness profiles for each galaxy. 
Average surface brightness profiles per morphological bin were created in order to 
characterize the variety of galaxy light profiles across the Hubble sequence. The 
ID bulge-disk decomposition parameters, as well as non-parametric galaxy mea- 
sures, such as effective radius, effective surface brightness and light concentration, 
are presented for all 286 VCC galaxies in each of the five optical/near- infrared 
wavebands. The profile decompositions account for bulge and disk components, 
spiral arms, nucleus and atmospheric blurring. The Virgo spiral galaxy bulges 
typically have a Sersic index n ~ 1, while elliptical galaxies prefer n ~ 2. No 
galaxy spheroid requires n > 3. The light profiles for 70% of the Virgo elliptical 
galaxies reveal the presence of both a spheroid and disk component. A more 
in-depth discussion of the structural parameter trends can be found in McDon- 
ald et al. (2009b). The data provided here should serve as a base for studies of 
galaxy structure and stellar populations in the cluster environment. The galaxy 
light profiles and bulge-disk decomposition results are available at the Centre de 
Donnees astronomiques de Strasbourg (CDS; http://cds.u-strasbg.fr/) and 
the author's own website (http://www.astro.queensu.ca/virgo). 



1. Introduction 

The advent of large-area near-infrared (NIR) surveys such as the Two-Micron All Sky 
Survey (Skrutskie et al. 2006; hereafter 2MASS) and the UKIRT Infrared Deep Sky Survey 
(Lawrence et al. 2007; hereafter UKIDSS) has forever improved our understanding of galaxy 
structure and evolution. The NIR data provide an unobscured view of galaxy structure, due 
to the relative insensitivity of NIR light to dust extinction, while probing the oldest and 
most representative stellar populations by mass. Unfortunately, the aforementioned surveys, 
designed largely for the study of bright infrared sources, use short exposure times (7.8 and 40 
seconds for 2MASS and UKIDSS, respectively), and suffer from the tremendous brightness of 
the NIR sky (typically ~ 3 orders of magnitude brighter than at r-band), leading to rather 
shallow surface brightness limits for galaxy studies when compared to similarly-designed 
optical surveys. 

The stability of the optical night sky, especially at r- and z-bands, allows for deep, 
accurate optical surface photometry with relatively short integrations and careful data pro- 
cessing. The Sloan Digital Sky Survey (Adelman-McCarthy et al. 2008; hereafter SDSS) has 
revolutionized optical astronomy, increasing the wealth of available astronomical data (both 
photometric and spectroscopic) by orders of magnitude. A development made possible by 
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the SDSS was the unraveling of bimodalities in various galaxy properties such as color, star 
formation rate and clustering (e.g., Strateva et al. 2003; Blanton et al. 2003; Kauffmann 
et al. 2003; Baldry et al. 2004; Brinchmann et al. 2004; Balogh et al. 2004). These and other 
studies have provided firm evidence that galaxies can be separated into two distinct classes: 
blue, star-forming galaxies found primarily in the field, and red, quiescent galaxies found 
primarily in clusters. 

Prior to these findings, Tully & Verheijen (1997; hereafter TV97) reported a different 
form of galaxy structural bimodality. This one, independent of the SDSS colour bimodality, 
was based on a sample of 62 members of the Ursa Major (UMa) cluster of galaxies; TV97 
found a bimodal distribution in the distribution of K'-band (2.2 /xm) surface brightness of 
galaxy disks. McDonald et al. (2009a), who re-analyzed the same UMa data using more 
sophisticated light profile decompositions to measure the disk structural parameters, corrob- 
orated TV97's results. Bell & de Blok (2000) have however argued that the UMa sample 
was too small for a statistically significant proof of bimodality. It is partly this unresolved 
issue, as well as the need for a complete NIR survey of a nearby galaxy cluster to study the 
(dust-free) structural parameters and stellar populations of cluster galaxies, that motivated 
our NIR survey of the Virgo cluster. 

Using the optical and NIR data presented in this paper, McDonald et al. (2009b) per- 
formed ID bulge-disk decompositions for 286 Virgo cluster galaxy surface brightness profiles 
and bolstered the evidence for a disk surface brightness bimodality. Furthermore, it was 
found for galaxies of all morphologies, that the effective surface brightness (defined as the 
surface brightness measured at the half-light radius) distribution had three peaks, each 
defining: (i) high surface brightness, gas-poor galaxies, (ii) low surface brightness, gas-poor 
galaxies and high surface brightness, gas-rich galaxies, and (iii) low surface brightness, gas- 
rich galaxies. A more thorough discussion of these, and other results from this survey, can 
be found in McDonald et al. (2009b). 

This paper presents the deep, optical and NIR surface brightness profiles for our sample 
of 286 Virgo cluster galaxies. The paper outhne is as follows: we describe in §2 the sample 
selection and the optical and NIR data acquisition. In §3, we give an overview of our NIR data 
reduction methods and we describe in §4 the process of surface brightness profile extraction 
for both the optical and NIR data. We discuss in §5 the quality of the NIR data and compare 
our data products with those from 2MASS and SDSS. The full collection of surface brightness 
profiles is presented in §6, along with a brief analysis of the average profile shape for various 
galaxy types. Galaxy structural parameters from our light profile decompositions and from 
non-parametric galaxy measurements are presented in §7. We conclude and ponder future 
work in §8. 
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In this paper, we assume a distance of 16.5 Mpc or m-M=31.18 for all Virgo cluster 
galaxies (Mei et al. 2007). At that distance, I" = 80 pc. 

2. Virgo Sample 

The construction of an unbiased distribution of galaxy surface brightness requires vol- 
ume completeness, which can be most easily obtained for galaxies in a cluster which all lie 
at a common distance. Our sample is drawn from the Virgo cluster catalogue of Binggeli 
et al. (1985; hereafter VCC). The VCC catalogue contains 2096 galaxies within an area of 
~140 deg2 on the sky, centered on the galaxy M87 at a~12'^25"' and (5~13°. The VCC 
is asserted to be complete down to a limiting absolute magnitude of Mb ~ — 13 and to 
contain many objects as faint as Mb ~ —11. The full sample was reduced to a manageable 
size of 286 galaxies by performing a brightness cut {Mb < —15.15 mag), a spatial cut (to 
remove contamination from W, W, and M background groups) and a velocity cut (K-ad < 
3000 kms~^). These cuts, which yield a sample covering a wide range of luminosities and 
morphologies, are discussed in further detail in McDonald et al. (2009b). 

2.1. NIR Data 

Deep iJ-band imaging for some VCC galaxies is already available from the 2MASS and 
GOLDMine0 (Gavazzi et al. 2003) databases. H-hand images from GOLDMine were kindly 
provided by G. Gavazzi, while calibrated 2MASS galaxy images were extracted from the 
online database. Many of the 2MASS and GOLDMine images were not deep enough for 
our purposes. Whilst adequate for large HSB galaxies, the high 2MASS surface brightness 
threshold (typically fin = 21 mag arcsec"^ - Bell et al. 2003; Courteau et al. 2007; Kirby 
et al. 2008) undermines the use of that database for deep extragalactic studies. Likewise, 
only a fraction of the available GOLDMine images were deep enough to properly separate 
bulge and disk light components, or suffered from ghost images. Whenever possible, we 
attempted to salvage such images by either masking the ghosts or confining our analysis to 
a single chip, which were typically unaffected by contamination due to ghosts. We secured 
if -band imaging for the remainder (187 galaxies) of our sample with the detectors ULBCam 
(UH 2.2-m), WFCAM (UKIRT), and WIRCAM (CFHT) over the period 2005-2009. The 
existing and new observations are summarized in Table 1 of McDonald et al. (2009b). 



http : / / goldmine . mib . inf n . it/ 
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The surface brightness profiles for the deep GOLDMine and 2MASS images were mea- 
sured using the same techniques as our new NIR images to ensure uniformity for the entire 
database. Further details regarding the data reduction process and quality are given below. 

2.2. Optical Data 

We have extracted calibrated ugriz images from the SDSS for the 286 VCC galaxies in 
our sample. Surface brightness profiles and total luminosities were extracted for all of these 
galaxies in all five SDSS bands using procedures described in Courteau (1996) and McDonald 
et al. (2009b). w-band images were consistently shallower than those in the griz bands and, 
thus, were discarded. Sky levels for background subtraction and the photometric zero- 
points for calibration were obtained from the SDSS image headers and the SDSS archives, 
respectively. However, we will see in §4.1| that we ultimately prefer our own estimates of 
the sky background for each of the galaxy images. The remainder of the profile extraction 
technique is identical to that used for the NIR photometry, as we describe below. 

3. NIR Observations and Data Reduction 

Virgo cluster galaxies which lacked suitable NIR imaging were observed with ULB- 
CanQon the 2.2-m UH telescope (Hall et al. 2004), WIRCAM on the 3.6-m CFHT (Puget 
et al. 2004) and WFCAM on the 3.8-m UKIRT (Hambly et al. 2008) from April 2005 to May 
2008 (see Table 1 of McDonald et al. 2009b). 

For all three NIR imagers, we only used one of the available detectors given the relatively 
small size of our targets. A single ULBCam detector has a 8f5 x 8^5 field of view (FOV) 
with a pixel scale of O'.'25/pixel, while WFCAM has a 13^6 x 13!6 FOV with a pixel scale of 
O'.'40/pixel, and WIRCAM has a 10^2 x 10^2 FOV with a pixel scale of O'.'30/pixel. Altogether, 
129, 35, and 27 galaxies were observed at H-band with ULBCam, WIRCAM and WFCAM, 
respectively. 

Given the rapid NIR sky fluctuations, we used a maximum single frame exposure time 
of 40s with ULBCam, 15s with WIRCAM and 10s with WFCAM to maximize sky flux, 
while keeping within the detector's linear regime. A dithering script was used to minimize 
resampling of bad pixels. In order to ensure sampling of each galaxy out to several effective 



•^ULBCam, the Ultra Low Background Camera, is a JWST prototype camera developed at the IfA by 
Don Hall and his team. It combines 4 2048x2048 arrays, of which we used only the cleanest one. 
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radii, we first pre-classified the galaxies by eye as high, medium, and low surface brightness. 
For these three classes we aimed for total exposure times of 8, 16, and 24 minutes (on 
the 2.2-m), respectively. These different exposure times were motivated by our desire for a 
well-sampled disk in order to obtain a reliable bulge-disk decomposition. 

Basic fiat-fielding, stacking and bad pixel rejection procedures were applied to the UL- 
BCam data using the XVista astronomical software packag^ Primary data reductions for 
the WIRCAM and WFCAM made use of observatory pipelines. The removal of geometric 
distortion and photometric calibration of the ULBCam data, as well as the procedures for 
extraction of surface brightness profiles for all NIR and optical data, are described below. 

3.1. ULBCam Geometric Distortion 

We corrected the significant geometric distortions in ULBCam images by observing the 
dense star field FS17 and computing the offset between the predicted and observed stellar 
positions, following Meurer et al. (2002). The following simple correction minimizes the 
distortion in our images: 

Ax{pix) = 16.64 - 0.0162?/ 
Ay{pix) = 0.0182X - 0.0230y + 3.80 

Additionally, each ULBCam 2048x2048 detector is comprised of four smaller 512x2048 
arrays that required horizontal offsets of -4 pixels (arrays 1 and 3) to account for their 
physical separation prior to determining the detector-wide astrometric correction. The 2D 
pattern of this correction model, shown in Fig. [l| is accurate to within 1 pixel (0.25") across 
the full FOV. 



3.2. NIR Flux Calibration 

The target flux calibration is based on 2MASS foreground stars. For each galaxy field, 
all stars above a given flux level (corresponding to the minimum detection level in 2MASS) 
were marked and their coordinates cross-correlated with the 2MASS stellar library. Typi- 
cally, three to seven 2MASS stars per ULBCam galaxy field were marked. The ratio of our 
instrumental magnitudes for those stars to those provided by 2MASS enabled a calibration 
of all of our photometry at each pixel, independent of any airmass or photometric variations 



" http : // astronomy . nmsu . edu/holtz/ xvista/ index . html 
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(our NIR images were still all obtained in most favorable conditions). This approach was 
tested extensively against the traditional method of multiple standard star observations per 
night and found to be equally rehable provided that at least three 2MASS stars were used 
per field. Under such conditions, the standard deviation of the derived zero-point corrections 
from each star is ~0.1 mag, as seen in Fig. [2j This error is comparable to those quoted for 
NIR photometric calibrations using standard star fields (de Jong 1996; Gavazzi et al. 1996; 
MacArthur et al. 2003). 

4. Surface Brightness Profile Extraction 

The measurements of galaxy surface brightnesses rely on the careful mapping, in one 
or two dimensions, of a representative light profile. Our profile extraction uses numerous 
XVISTA routines and procedures described in Courteau (1996). We review the important 
steps below. While azimuthally-averaged surface brightness profiles for some of our galaxies 
are available elsewhere (e.g. 2MASS and GOLDMine), we have recovered original images in 
all cases and rederived our own surface brightness profiles for all of these galaxies, to ensure 
a homogeneous database. 

4.1. Sky Measurement 

A sky level error of only ±0.01% at NIR wavelengths can result in either an artificial 
truncation (over-subtraction) or upturn (under-subtraction) of the galaxy SB profile at large 
radii. At H-band, this truncation/upturn would occur at ~ 23.5 mag arcsec"^. An accurate 
assessment of the sky level is thus crucial for the extraction of a reliable, deep, NIR surface 
brightness profile. 

For all images, the sky level is estimated in the final fully-reduced image by isolating five 
regions (by eye) away from the galaxy, free of any other sources, and calculating the mode 
of the sky intensities per pixel within each sky region. The average and standard deviation 
of the five sky values were then computed. Sky levels were estimated this way not only for 
our data but for 2MASS and GOLDMine images as well. The typical sky fluctuation for our 
NIR data is -0.005%. 

Thanks to its extended sky coverage, SDSS should achieve accurate sky level measure- 
ments (http : //www. sdss . org/dr5/algorithms/f lat_field.html). However, we still use 
our interactive assessment of the sky background, as described above, for all SDSS data. 
Fig. |3] shows the difference between our measured sky values and those provided by SDSS. 
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The latter is always biased high, likely due to the inclusion of bright, extended sources. Our 
interactive technique ensures a (mostly) contaminant-free selection of sky fields. However, 
due to the fact that we extend the surface brightness profiles to the edge of the field of view, 
the uncertainty in the profile shape is dominated by systematic (sky measurement) errors 
at large radii. While seemingly small, sky level errors of a few percent, as seen in Fig. [3} 
yield significant deviations in a galaxy surface brightness profile (McDonald et al. 2009b). 
We find that the shape of the surface brightness profile in g-hand is typically dominated by 
sky measurement errors by fig ~26 mag arcsec"^. Below this surface brightness level (and 
corresponding surface brightnesses in r,i,z) users should consider the data with 1-sigma sky 
error envelopes. 

4.2. PSF Measurement and Star Masking 

The seeing point-spread functions (PSF) is often modeled either as a Gaussian or a 
Moffat (1969) function. The latter has broader wings, which better matches our stellar fields; 
we have thus used the Moffat function throughout our analysis. On average, 17 stars per field 
were modeled for PSF measurements. The average FWHM for all ULBCam observations 
taken over 2005-2008 was 1.2 ± 0.2"; slightly better imaging was achieved at CFHT and 
UKIRT. The average FWHM measurement for each image is later used to convolve the 
analytical models for the galaxy bulge and disk light (see McDonald et al. 2009a). 

The identification of foreground stars also enables their removal from the galaxy light. 
The use of circular masks with radii equal to 4x FWHM ensures that most of this contam- 
inant light is properly removed. Our automated masking routine could not identify bright 
stars with diffraction spikes or other irregularly shaped features; those obtrusive objects were 
manually masked. 

4.3. Isophotal Fitting 

Surface brightness profiles were extracted for all of our galaxy images. The XVISTA 
command, PROFILE, performs this operation through isophotal fitting, using a generalized 
non-linear least-squares fitting routine. The ellipticity and position angle, but not the center, 
of each elliptical isophote are allowed to vary. The center is determined by measuring the 
centroid of the central light distribution with the XVISTA command AXES. Beyond 1-2 
disk scale lengths, the position angle and ellipticity usually settle to constant values. The 
isophotal solution based on these values is extrapolated to larger radii where the signal-to- 
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noise is too low for PROFILE to converge on a unique solution. Adjustments to the ellipticity 
profile can be made with the XVISTA command, CPROF, to account for abrupt changes due 
to any non-axisymmetric features. Ellipse fitting to map the radial surface brightness profile 
requires the galaxy to be slightly (but not fully!) inclined, so that the azimuthal direction is 
projected onto the plane of the sky. Edge-on galaxies are thus excluded from this procedure. 
Fortunately, only 3 galaxies in our sample have inclinations > 80°; non-parametric structural 
measures (e.g., effective surface brightness) for those galaxies are still valid. 



4.4. Profile Depth and Signal-to-Noise Ratios 

We calculate the signal-to-noise (hereafter S/N) ratio of a surface brightness profile as 
a function of radius as: 

where the number of pixels (expressed as an area) along each isophote is given by: 

= 27rv/0.5(a2 + 62), (2) 

where a, b are the semi-major and semi-minor axes of a given elliptical isophote, It{r) is its 
total surface brightness level (intensity in counts), and Isky is the surface brightness of the 
sky. 



The measurement of Igky was outlined in ^4.1 According to Eq. ([T]), depths of ~ 23.5 H 



mag arcsec"^ with S/N ~3 can be achieved with an exposure time of 480s (at the 2.2-m 
UH telescope). For intermediate and low surface brightness galaxies, our exposure times of 
960s and 1440s were determined to yield S/N 3 at surface brightness depths of 24 and 
24.5 H-mag arcsec"^ , respectively. These levels correspond to ~4 disk scale lengths for disk 
galaxies and ~5 effective radii for spheroidal systems. 



5. NIR Data Quality 

While our collection of Virgo cluster galaxy light profiles results from the merging of 
independent imaging surveys, we have imposed our own uniform analysis methods to the 
entire database. Internal errors can be assessed from repeat measurements, and we can 
compare our extracted light profiles with those from the original authors, whenever available. 
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We can also compare our estimates of empirical quantities such as total luminosity, scale radii 
and concentration to those provided with surveys such as 2MASS and SDSS in order to test 
for systematic differences. 

5.1. Independent Calibration Errors from Multiple Measurements 

We have repeat measurements for three galaxies observed from 2005 to 2007: VCC1614, 
VCC1516 and VCC1664. We find no variation in the surface brightness profiles over these 
three epochs, demonstrating the stability of our calibrations over time. Since our calibration 
of ULBCam data is tied to the 2MASS system, we did not derive extinction coefficients and 
zero-point corrections for each observing night. However, we can compute relative offsets 
between the observed and 2MASS magnitudes for several stars in each observed field. The 
distribution of the photometric offsets shown in Fig. |4] spans 4 years of observations and a 
full range of airmasses. The measured width of ~ 0.2H mag for this distribution is expected 
for night-to-night zero-point variations alone (e.g., Courteau 1996). 

5.2. Profile Comparison with Independent Measurements 

A total of 80 GOLDMine galaxies that satisfy our Q criterion (McDonald et al. 2009b) 
were included in our sample. Original images from GOLDMine's heterogeneous collection of 
galaxies were kindly provided to us by G. Gavazzi. The average if-band surface brightness 
zero-point offset of ± 0.15 mag arcsec"^ between our ULBCam and GOLDMine's original 
light profiles is small enough to justify merging the two samples into one. However, for 
complete uniformity, we have recomputed surface brightness profiles for GOLDMine galaxies 
using our own data reduction techniques outlined earlier. Another motivation for doing so is 
the desire to have linearly, rather than logarithmically (as in GOLDMine), sampled surface 
brightness profiles. This exercise also led to more data points in the outer disks, which is 
relevant for the extrapolation of light profiles to infinity. Fig. |5] shows the overlap of our 
profiles derived from ULBCam and GOLDMine images for selection of VCC galaxies covering 
a broad range in luminosity and morphology. 

We now compare our calibrated profiles with the fully homogeneous database of deep H- 
band images of late-type northern spirals by MacArthur et al. (2003; hereafter M03). These 
if-band images of spiral galaxies were collected at the KPNO 4-m telescope with the COB 
detector; typical total integration times were 20 mins per galaxy. We observed 12 non- Virgo 
UGC galaxies from MOS's sample with ULBCam in April 2007 for comparison. As shown in 
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Fig. [6| the match with M03 is typically very good, with a mean deviation between the two 
samples of ~ 0.1 mag arcsec"^. This comparison with an independent database reinforces 
the quality of our dataset. 

To summarize this section, we have shown that we can merge our Q > 0.5 GOLD- 
Mine and ULBCam if-band SB profiles into one sample with a systematic error < 0.15 
mag arcsec"^. Because all profiles are calibrated to 2MASS, the 2MASS profiles with Q > 0.5 
are also naturally integrated into our system. Our entire collection of iJ-band SB profiles is 
thus self-consistent to within < 0.15 mag arcsec"^ out to the last measured data point. 

5.3. Comparison of 2MASS and SDSS Data Products 

We can also compare our well-tested, non-parametric data products, such as the half- 
light radius and total magnitude, with those provided by automated 2MASS and SDSS 
pipelines. The upper panel of Fig. [7] shows our total, extrapolated H-hand brightnesses 
versus those provided by 2MASS. The agreement between the two samples for H < 12 mag 
is excellent, however the shallower 2MASS catalog is biased for fainter magnitudes. The 
lower panels of Fig. [7] show our H-hand radii against the 2MASS i^'-band isophotal radius, 
rK2o, measured at the i^-band surface brightness level of 20 mag arcsec"^. We compare rK2o 
to two radii, r^o and rgo, which enclose 50% and 80% of the total H-hand light, respectively. 
Besides the bandpass differences, the scatter in this distribution is largely due to the bias 
between an absolute metric, rx2o, and a relative metric, r^Q or rgo; the former is fixed for any 
galaxy profile whereas the latter shifts radially as a function of galaxy profile shape. The 
two measures thus differ as a function of galaxy mean surface brightness, as shown in Fig. [7j 
the agreement worsens for fin > 18 mag arcsec"^ . Overall, for tuh < 13 mag and fin < 18 
mag arcsec"^ , the 2MASS total magnitudes and isophotal radii are reliable to within ~ 0.1 
mag and 0.2 dex, respectively. 

Before performing a similar comparison with the SDSS data, we must first confirm that 
we have properly extracted and calibrated these data. To that end, we extract Petrosian 
magnitudes from our cahbrated profiles and compare these to the published SDSS Petrosian 
magnitudes (Adelman-McCarthy et al. 2008). For a fair comparison, we use the distributed 
Petrosian radii and computed the total brightness within a circular aperture of that radius. 
The results from this analysis, shown in Fig. |8} confirm that we reproduce the published 
Petrosian magnitudes to within their error for nearly all galaxies. This confirms our proper 
calibration of the SDSS surface brightness profiles. 

A comparison with various SDSS data products is shown in Fig. [9j The top figure 
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shows the SDSS Petrosian r-band magnitude versus our total isophotal r-band magnitude 
from SDSS images. Many galaxies have Petrosian magnitudes more than 2 magnitudes 
fainter than ours, likely due to a galaxy misidentification by the SDSS pipeline (see Hall 
et al. 2011, in prep.). The mean ~ 0.2 mag offset from the 1-to-l line for well-behaved 
systems is simply the reflection of comparing our total magnitudes extrapolated to r = oo 
with Petrosian magnitudes which are measured within a finite aperture. The deviation at the 
bright end is due to the SDSS pipeline "shredding" of large galaxies (e.g., Panter et al. 2004). 
The large scatter at all magnitudes is most likely due to the misidentification of clumpy, star- 
forming galaxies into several individual sources. The bottom left of Fig. |9] shows the SDSS 
Petrosian half-light radius, rs^sDss, against the half-light radius, rso, measured with our own 
software from SDSS images. The deviation from the 1-to-l line is largely due to comparing 
radii derived from circular (SDSS) apertures versus elliptical apertures, since, at a given 
radius, a circular aperture will always enclose more flux than an elliptical aperture of the 
same maximal size. 

Finally, the lower-right corner of Fig. [9] shows a comparison of our and the SDSS con- 
centration parameter, defined as Cxy=51ogio(ry/rx), where rx and vy contain X x 10% 
and Y x 10% of the total light of the galaxy. The curvature between the radii seen in the 
lower-left panel of Fig. [9] cancels out to yield a rather linear, albeit noisy, mapping between 
the two concentrations. The conversion between the two concentration parameters is given 

by C59 = O.34C28 + 0.71. 

Overall, these comparisons highlight the risks of taking survey data at face value (see 
also Hall et al. 2010). 



6. Surface Brightness Profiles 



We show in Fig. 10 the griz and i7-band surface brightness profiles for 12 of the 286 



VCC galaxies in our main sample (see http://www.astro.queensu.ca/virgo/ for the full 



collection). The full morphological range spanned by our sample from blue to red and 
dwarf to giant is demonstrated by the wide variety of light profile shapes. To demonstrate 
this, we have computed the average r-band surface brightness profile for galaxies in several 



morphological bins. Fig. 11 shows the differences and similarities between light profiles at 
all morphological types, normalized at yUg and rg. In a plot of /i — /ig versus r/rg, all SB 
profiles must pass through the point (1,0). In each window we show the normalised surface 
brightness profiles of a given morphology, with the mean profile per morphological bin shown 
in color. The middle and lower windows to the right show the average surface brightness 
profiles in each morphological bin interior and exterior to r^,. For r < O.Srg, gas-poor (E, 
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SO, Sa) galaxies exhibit classical cuspy bulges, while gas-rich types (Sb, Sc, Sd, dwarf and 
irregular galaxies) have cored profiles, commonly associated with pseudo-bulges (MacArthur 
et al. 2003; Kormendy & Kennicutt 2004). For r > 0.3re, the relative shape of the outer 
profiles can be compared. Both giant and dwarf spheroidal systems (E, dE, SO, dS) have 
similar shapes, while disk-like systems (Sa-Sd) exhibit a more exponential outer profile. 

Similar conclusions are reached if we bin the profiles by concentration rather than mor- 



phological type (Fig. 12). The use of a quantitive measure for the concentration yields 
tighter distributions as opposed to the more subjective morphological classification. There 
is a natural division for systems with pseudo bulges (C<3.8) versus those with a classical 
bulge (C>3.8), as shown by McDonald et al. (2009b). 

We find the same results with the H-band data, though the r-band data are, typically, 
deeper than at H-band. The most noticeable difference is with the depth of the dwarf and 
irregular galaxy profiles. The average dE, dS and Irr profiles are well defined out to much 
larger radii at r-band than at if-band. Overall, the trends that are present in the mean 
i^-band profiles are reproduced at r-band, suggesting that these mean profiles are robust. 

Whereas Lauer et al. (2007) hinted at a bimodal separation of the central regions of 
giant and dwarf ellipticals into cuspy and cored profiles, respectively. Cote et al. (2007) 
showed that the distribution of central profile shapes is a continuum. Cote et al. suggested 
that a possible selection bias and choice of radius at which the inner slope is measured could 
explain the differences, showing that a continuum of inner slopes is seen in the comparison of 



the inner profile slope with a non-parametric parameter such as total luminosity. Our Fig. 12 
attests to the continuous distribution of inner slopes when we bin by galaxy concentration 
(which correlates with total luminosity). Thus, a dichotomy between cuspy and cored profiles 
is only seen if profiles are binned by morphology; a continuous distribution prevails with a 
more physically motivated binning criterion such as concentration. 



7. Galaxy Structural Parameters 

In order to determine the contribution of the bulge and disk components to the total 
galaxy light, we have performed 1-D bulge-disk decompositions on the full sample of Virgo 
cluster galaxies following the techniques described in M09a. Due to the 1-D nature of these 
decompositions, we are unable to fully account for non-axisymmetric components such as bars 
and ovals. These components will necessarily be added to either the bulge or the disk light, 
adding a degree of uncertainty to these fits. Ideally, one should use a 2-D decomposition 
code such as GIM2D (Simard et al. 2002), BUDDA (de Souza et al. 2004) or GALFIT 



- 14 - 



(Peng et al. 2010) so that non-axisymmetric features can be properly modeled. For further 
discussion of the relative merits of 1-D and 2-D bulge-disk decompositions, see MacArthur 
et al. (2003). While 2D decompositions can reduce some of the inherent uncertainty in 
ID model parameters, differences between ID and 2D values are at the 10-20% level. We 
proceed with our ID approach here; our 2D decompositions will be presented elsewhere. We 
describe the salient features of this method below. 

Our bulge model uses a generalized Sersic function (Sersic 1968) with three free param- 
eters: an effective (or half-light) radius, rg, the surface brightness at that radius, Jg, and a 
shape parameter, n. For the special case of n = 1, the Sersic function reduces to a simple 
exponential. The disk light is modeled with an exponential function described by two free 
parameters: a disk central surface brightness, /q, and a disk scale length, h. The bulge and 
disk models are convolved with the Moffat (1969) function to simulate the effects of seeing 
on the intrinsic light profile of a galaxy. 

For spiral and irregular galaxies (SO-Sd, Irr) we model the surface brightness profile 
with a combination of a Sersic bulge and exponential disk. If the inclusion of the bulge 
component does not sufficiently improve the fit quality (i.e. Sd or Irr galaxies), the software 
will discard it and fit an exponential function alone in an attempt to decrease the degrees of 
freedom and, thusly, minimize the reduced x^- Since spheroidal galaxies (E, dE, dS) may or 
may not have a disk, we also attempt to model the profile with a single Sersic. Once again, if 
the absence of a disk component increases the goodness-of-fit, the software will automatically 
discard this component of the model. The fit minimization uses brightnesses in magnitude 
(not counts), and the lowest global value of the per degree of freedom determines the best 
fit. 

Our model decomposition can also account for the presence of a nucleus, here modeled 
by a seeing-convolved delta function at r = (one additional free parameter, rrinuc), as well 
as spiral arms, which are modeled by a smooth increase in brightness above the underlying 
disk. We find that the addition of a nuclear component reduces the need for high-n bulges, 
yielding an average shape parameter of n ~ 1. These additional features increase the degrees 
of freedom and, thus, are only applied for obvious cases. 

Non-parametric structural measurements, such as the total magnitude (extrapolated 
from the last data point to infinity), the effective radius and surface brightness, and galaxy 
fight concentration, were computed for all 286 galaxies in each of the 5 bands. These mea- 
surements are nearly independent of the bulge-disk model decompositions, relying only on 
the fit in order to determine the amount of galaxy light beyond the last data point (typically 
< 0.2 mag). The fact that these parameters arc independent of any model (i.e. no assump- 
tion that a galaxy profile can be separated into components) makes them most useful for 
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galaxy structure studies. The effective radius, rg, is defined as the radius which encloses half 
the total light, while the effective surface brightness is the surface brightness at this radius. 



He = fJ'ire). We compute the concentration parameter, C28, as in [5.3 



The parametric and non-parametric structural measurements for our 286 VCC galaxies 
in the griz and H bands are provided online at the CDS ( http : //cds . u-strasbg . f r/[ ) and 
our private website (http : / / www . astro . queensu . ca/virgo/) . For spiral galaxies with well- 
defined bulges and disks (i.e., Sa-Sb), we find typical r-band errors on our fitting parameters 
of: A/ie,d = 0.1 mag arcsec"^ , Afie,b = 0.2 mag arcsec"^ , Ave^d = 3%, Are,fe = 15%, and 
Arib = 0.3. These uncertainties are similar to those reported by MacArthur et al. (2003), who 
relied on similar techniques. For irregular or poorly-sampled systems, these uncertainties can 
increase by as much as an order of magnitude. The distribution of the H-band parametric 
quantities is shown in Fig. 13 for different morphological bins. This figure clearly shows 
the disk surface brightness bimodality discovered by TV97 and confirmed by McDonald 
et al. (2009a,b). These distributions are meant to be compared to galaxy formation models. 



8. Tables of Median Structural Parameters 

We provide in Table 1 the median and standard deviation of the four non-parametric 
parameters (total magnitude, effective surface brightness, effective radius, and concentra- 
tion) as a function of morphology and bandpass. This table shows the variation of galaxy 
structural parameters with wavelength and morphology, as well as the spread of each pa- 
rameter within a given morphological/wavelength bin. For instance, there is a substantial 
range in elliptical galaxy brightnesses, as expected from the bimodal distribution of ellipti- 
cal surface brightnesses (McDonald et al. 2009b), and a clear division of galaxies into cuspy 
(C28 ~ 4.0) and disk-like (C28 ~ 2.8) galaxies, also as reported by McDonald et al. (2009b). 
Additionally, Table 1 shows that disk galaxies have, on average, effective radii that are twice 
as large as spheroidal galaxies, regardless of surface brightness or color. 

Table 2 gives the median and standard deviation of five parametric quantities from ID 
bulge-disk decompositions. This includes three bulge parameters: fj,e,b, re,b, n and the two 
disk parameters: /ie,d and r^^d- This table shows the variation of bulge and disk structural 
parameters as a function of color and morphology. We find that the "typical" bulge shape 
is indeed exponential (e.g., Courteau et al. 1996; MacArthur et al. 2003) and that the light 
profile of a typical elliptical galaxy is well fit by a single Sersic index of n ~ 2. No Virgo 
galaxy spheroid requires n > 3. Futhermore, the decomposition of elliptical galaxy light 
profiles requires the addition of a disk component for ~ 70% of our sample. These disks 
tend to have low surface brightness, with average brightnesses similar to Sc-Sd galaxies, but 
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typically have scale lengths twice as large. Finally, we find that the typical spiral galaxy 
disk has a half-light radius ~ 5x larger than that of the bulge. 

Tables 1 and 2 offer a wealth of information about the structural parameters of galaxies 
across the Hubble sequence and over a range of wavelengths. An extensive discussion of 
these trends is found in McDonald et al. (2009b). These tables ought to provide a valuable 
compendium of galaxy structural properties for theoretical investigations of galaxies as well 
as a much needed foil for comparisons with high redshift galaxy samples. 



9. Summary 

We have presented the results of a NIR survey, supplemented with optical imaging from 
the SDSS, of 286 morphologically-diverse Virgo cluster galaxies. Our NIR images come from 
various sources including archival if-band data from the 2MASS and GOLDMine archives 
and new H-hand data for 171 galaxies from the UH 2.2-m, CFHT and UKIRT telescopes. 
These NIR data have been carefully calibrated, yielding a photometric accuracy of ~ 0.1 
mag arcsec"^ over the entire sample. Both the new and archival data were analyzed in a 
homogeneous manner using our own software for uniformity. For each galaxy in our sam- 
ple we have extracted g,r,i,z,H surface brightness profiles yielding a deep, multiwavelength 
collection of surface photometry. Full bulge-disk decompositions were performed in each 
of the five wavebands for all 286 VCC galaxies. We also provide the average structural 
parameters for galaxies in various morphological bins, offering a most valuable benchmark 
for comparison against theoretical investigations of galaxies as well as high redshift galaxy 
samples. This database has recently been utilized in studies of galaxy structural param- 
eters by McDonald et al. (2009b) and stellar populations by Roediger et al. (2010) and 
Prochaska et al. (2011). The galaxy light profiles and results of bulge-disk decomposi- 



tions are available online at the CDS (http://cds.u-strasbg.fr/) and our own website 
(http : / / www . astro . queensu . ca/virgo/) . 
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Fig. 1. — Geometric distortion in one array of the ULBCam 4096x4096 detector. The 
magnitudes of the vectors have been increased by a factor of 8 for clarity. 
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Fig. 2. — Calibration error as a function of the number of field stars used per calibration. 
In order to compute the standard deviation, N stars were chosen at random and the mean 
brightness was computed. This process was repeated 10,000 times to ensure every possible 
combination of stars was used. The standard deviation represents the difference in means 
depending on which N out of 11 stars were chosen. 
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Fig. 3. — Difference between our interactiveiy-measured and tlie SDSS-provided sky fiux as 
a function of our measured sky level (in counts). Negative Y- values indicate a larger SDSS 
value. 
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Number of Stars 

Fig. 4. — Upper panel: Distribution of magnitude offsets between observed and 2MASS 
measurements for foreground stars over many observing runs. The black, red and blue 
histograms correspond to observations taken in 2005, 2006 and 2007, respectively. Lower 
panel: Number of foreground stars used to determine a relative magnitude offset per single 
observation; a minimum of 3 stars per field was found to be necessary to constrain the 
photometric zeropoint. 
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Fig. 5. — Comparison of GOLDMine (red) and ULBCam (black) i^-band surface brightness 
profiles for a sample of matching galaxies, ordered (from left to right and top to bottom) by 
total luminosity. The good match in profile depth found at the bright end of the sample is 
lost for fainter galaxies. 
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Fig. 6. — Comparison of M03 (red) and ULBCam (black) i7-band surface brightness profiles. 
These galaxies, while not members of the Virgo cluster, were observed using ULBCam in 
order to provide an independent comparison of our calibration and analysis techniques. 



- 26 - 




Fig. 7. — (Top) Comparison of the 2MASS total H-hand magnitude with our total isophotal 
H-hand magnitude. Up to 12 if-mag, the agreement is superb with a typical 1-a deviation 
of 0.2 mag. The shallower 2MASS is biased towards fainter galaxy magnitudes for ttlh > 12 
mag. (Bottom) Comparison of the 2MASS i^'-band isophotal radius, rK2o, measured at the 
i^'-band surface brightness of 20 mag arcsec"^, with our own radii, Vfj^o and rj^go, which 
enclose 50% and 80% of the total iJ-band light, respectively. Modulo the different band- 
passes, the scatter is largely due to the different light profile shapes for galaxies with different 
mean surface brightness. The red points for < 16 mag arcsec"^ , the blue points for 
fin > 18 mag arcsec"^ , and green points for the galaxies in between, make this conclusion 
very clear. The scatter in the distributions of rK2o, rn^o or rnso, and surface brightness 
naturally increases for fainter galaxies. 
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Fig. 8. — Comparison of available Petrosian magnitudes from the SDSS pipeline and those 
calculated based on our surface brightness profiles extracted from SDSS images. The red, 
dashed lines represent the one-to-one relations. These plots confirm that our zeropoint 
calibration and magnitudes are properly computed. 
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Fig. 9. — (Top) Comparison of the SDSS Petrosian r-band magnitudes with our total isopho- 
tal r-band magnitudes from SDSS images. The offset at the bright end is due to the SDSS 
pipehne "shredding" large galaxies. The large scatter is also due to SDSS's misidentification 
of clumpy, star-forming galaxies into multiple sources. (Bottom) Left: Comparison of the 
SDSS Petrosian r-band half-light radii, rs^sDss, with the r-band half-light radii, rso, mea- 
sured with our own software from SDSS images. The deviation from the one-to-one line 
is largely due to comparing radii derived from circular (SDSS) apertures versus elliptical 
apertures. Right: Comparison of SDSS concentration, C59, with our own measure, C28- 
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Fig. 10. — Surface brightness profiles for a subset of the 286 galaxies in our sample. The 
colour coding is as follows: g (blue), r (green), i (orange), z (red), H (black). The VCC 
galaxy name is included in the upper right corner of each window, with the H-band data 
source beneath it. The data source codes are: (UH 2.2-m ULBCam observations from 

2005/06/07/08), CFH (CFHT WIRCAM), UKT (UKIRT WFCAM), 2MS (2MASS archive). 
All other codes (C##, N##, T##, G##) refer to data taken from the GOLDMine survey 
and correspond to the telescope used and year of the observation. (The SB profiles for all 



286 galaxies are available online at http://www.astro.queensu.ca/virgo/). 
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r— band Normalized Surface Brightness Profiles 




Fig. 11. — r-band surface brightness profiles for our entire sample, rescaled in terms of /ig 
and Te, and binned by morphological classes. Average profiles for each bin are represented by 
colored lines. The middle and lower windows in the right column show the average profiles 
interior and exterior to rg. 
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Fig. 12. — r-band surface brightness profiles for our entire sample, rescaled in terms of fie 
and Tg, and binned by concentration. Average profiles for each bin are represented by colored 
lines. The middle and lower windows in the right column show the average profiles interior 
and exterior to Vp. 
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Fig. 13. — Distribution of H-band structural parameters from bulge-disk decompositions for 
Virgo galaxies sorted into 6 different morphological bins. 
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Tabic 1: Median and standard deviation of Virgo cluster galaxy non-parametric structural 
quantities for different morphologies and bandpasses. 



Parameter 


Morpfiology 


9 


r 


i 


z 


H 




E 


13.5(1.8) 


12.9(1.8) 


12.4(1.8) 


12.7(1.9) 


10.6(2.2) 


[mag] 


SO 


12.5(1.3) 


11.6(1.4) 


11.2(1.4) 


11.0(1.5) 


9.3(1.6) 


(1) 


Sa-Sb 


12.1(1.2) 


11.5(1.3) 


11.1(1.3) 


10.9(1.3) 


9.0(1.5) 




Sc-Sd 


13.3(1.3) 


12.7(1.5) 


12.6(1.4) 


12.4(1.5) 


9.9(1.8) 




dE 


15.3(0.9) 


14.7(0.8) 


14.4(0.9) 


14.2(1.3) 


12.5(0.9) 




dS 


14 7(0 5"! 


Xt:. J. I U.U J 


1 3 R(() ft) 


1 3 fifn 7) 

XO.Ul U. / J 


1 1 7(0 q"! 




Irr 


14.9(1.3) 


14.3(1.4) 


14.0(1.5) 


13.7(1.5) 


12.3(1.6) 




E 


22.0(1.4) 


21.2(1.5) 


20.9(1.5) 


20.5(1.4) 


18.3(1.5) 


[mag arcsec"^] 


SO 


22.3(0.9) 


21.5(1.0) 


21.1(1.0) 


20.8(1.0) 


18.9(1.1) 


(2) 


Sa-Sb 


22.8(0.7) 


22.2(0.7) 


21.8(0.7) 


21.5(0.7) 


19.5(0.7) 




Sc-Sd 


23.5(0.7) 


23.0(0.7) 


22.7(0.7) 


22.3(0.7) 


20.6(1.0) 




dE 


24.1(0.7) 


23.4(0.7) 


23.1(0.8) 


22.6(0.7) 


21.1(0.9) 




dS 


23 R(0 Q) 


23 2(0 R) 


22 Qfn R) 

.U \\J j 


22 Ft(() R) 


20 fifn R) 

.\J\\J .\D J 




Irr 


23.7(1.0) 


23.3(1.1) 


23.1(1.1) 


22.6(1.1) 


21.3(1.3) 




E 


1.1(0.8) 


1.1(0.8) 


1.1(0.9) 


1.0(0.7) 


0.9(1.1) 


[kpc] 


SO 


1.7(0.9) 


1.8(0.9) 


1.7(0.9) 


1.6(0.8) 


1.4(1.0) 


(3) 


Sa-Sb 


2.6(1.6) 


2.5(1.5) 


2.9(1.4) 


2.5(1.4) 


2.5(2.0) 




Sc-Sd 


2.4(1.8) 


2.4(1.7) 


2.4(1.5) 


2.3(1.4) 


2.2(2.4) 




dE 


1.3(0.5) 


1.3(0.5) 


1.2(0.5) 


1.0(0.8) 


1.2(0.5) 




dS 


1.2(0.6) 


1.2(0.5) 


1.2(0.5) 


1.2(0.5) 


1.1(0.5) 




Irr 


1.4(0.7) 


1.4(0.7) 


1.5(0.8) 


1.3(0.7) 


1.3(0.8) 


C28 


E 


3.8(0.7) 


3.9(0.7) 


4.0(0.7) 


3.8(1.0) 


3.8(0.7) 


(4) 


SO 


4.0(0.8) 


4.1(0.8) 


4.1(0.8) 


4.2(0.8) 


4.0(0.9) 




Sa-Sb 


3.1(0.7) 


3.2(0.7) 


3.2(0.7) 


3.3(0.7) 


3.4(0.7) 




Sc-Sd 


2.6(0.4) 


2.5(0.5) 


2.8(0.5) 


2.7(0.6) 


2.8(0.4) 




dE 


2.9(0.6) 


3.0(0.6) 


2.9(0.6) 


2.7(0.5) 


2.9(0.5) 




dS 


3.3(0.4) 


3.3(0.4) 


3.2(0.4) 


2.9(0.4) 


3.3(0.6) 




Irr 


2.7(0.8) 


2.7(0.8) 


2.6(0.7) 


2.5(0.7) 


2.6(0.7) 



(1) Total apparent magnitude extrapolated to r=oo (mag) (2) Effective surface brightness (mag arcsec ^) 
(3) Effective radius (kpc) (4) Concentration: C28=51ogio(r8o/r2o) 
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Table 2: Median and standard deviation of Virgo cluster galaxy parametric structural pa- 
rameters from bulge-disk decompositions for different morphologies and bandpasses. 



Parameter 


Morphology 


9 


r 


i 


z 


H 


n 


El 


2.4(1.1) 


2.6(1.0) 


2.3(1.0) 


2.2(1.4) 


2.1(0.9) 


(1) 


E2 


1.4(0.6) 


1.5(0.5) 


1.5(0.6) 


1.3(0.7) 


1.5(0.7) 




SO 


1.4(0.6) 


1.4(0.5) 


1.4(0.6) 


1.3(0.6) 


1.2(0.7) 




oa-oD 


1 0{C\ '7 \ 
L.Zyj. ( j 


1 0(C\ f\\ 
i.Zt^U.D ) 


1 A(C\ (\\ 
i.4(^U.D ) 


1 0{C\ '7 \ 
i.ZI^U. ( ) 


1 fi/'n f{\ 

i.UI^U.O ) 




Sc-Sd 


0.5(0.3) 


0.7(0.5) 


0.6(0.4) 


0.6(0.4) 


0.8(0.5) 




dE/dS 


0.9(0.3) 


0.9(0.4) 


0.9(0.4) 


0.8(0.4) 


0.8(0.5) 




Irr 


0.6(0.6) 


0.8(0.6) 


0.7(0.5) 


0.5(0.5) 


0.7(0.6) 




El 


22.7(1.4) 


21.7(1.4) 


21.4(1.4) 


21.3(1.5) 


18.7(0.9) 


[mag arcsec"^] 


E2 


21.0(1.4) 


20.5(1.4) 


20.2(1.4) 


19.8(1.6) 


17.5(1.3) 


(2) 


SO 


20.6(1.3) 


20.0(1.3) 


19.4(1.4) 


19.2(1.5) 


17.2(1.3) 




Dcl-oD 










1 7 (\(^ ^\ 

i i .Oi^i.O j 




Sc-Sd 


22.6(1.5) 


21.9(0.9) 


21.6(1.0) 


21.4(1.2) 


19.7(1.4) 




dE/dS 


23.0(0.9) 


22.6(0.9) 


22.2(1.0) 


21.9(0.9) 


20.0(1.1) 




Irr 


22.6(1.4) 


22.1(1.5) 


21.8(1.6) 


21.0(1.7) 


20.2(1.8) 


re,b 


El 


1.7(2.3) 


1.6(2.0) 


1.5(2.0) 


1.5(3.2) 


1.3(2.0) 


[kpc] 


E2 


0.6(0.4) 


0.6(0.3) 


0.6(0.3) 


0.6(0.3) 


0.7(0.9) 




Qn 


u.ot^u.o ) 


u.o^^u.o ) 


U.Ot^U. ( ) 


U.DI^U.O ) 


u.oi^u.y J 




Sa-Sb 


0.5(1.4) 


0.6(0.6) 


0.6(0.8) 


0.6(0.7) 


0.5(0.6) 




Sc-Sd 


0.5(0.4) 


0.4(0.4) 


0.6(0.6) 


0.6(0.5) 


0.7(1.1) 




dE/dS 


0.6(0.4) 


0.6(0.3) 


0.6(0.3) 


0.5(0.3) 


0.5(0.4) 




Irr 


0.6(0.4) 


0.8(0.5) 


0.8(0.4) 


0.4(0.3) 


0.8(0.5) 




E2 


23.8(1.6) 


23.3(1.6) 


22.7(1.6) 


22.3(1.6) 


20.9(1.8) 


[mag arcsec"^] 


SO 


22.9(1.5) 


22.3(1.3) 


21.9(1.2) 


21.3(1.4) 


19.6(1.5) 


(4) 


Sa-Sb 


22.6(1.3) 


22.0(0.9) 


21.6(1.1) 


21.3(1.1) 


18.9(1.1) 




Sc-Sd 


22.9(0.8) 


22.4(0.9) 


22.1(1.1) 


21.9(0.8) 


20.5(1.4) 




dE/dS 


24.3(1.0) 


23.7(0.9) 


23.4(0.9) 


22.9(0.9) 


21.4(1.0) 




Irr 


23.8(1.2) 


23.4(1.1) 


23.2(1.1) 


22.6(0.9) 


21.4(1.3) 


re,d 


E2 


2.4(2.3) 


2.4(2.6) 


2.4(2.4) 


2.1(2.5) 


2.5(9.1) 


[kpc] 


SO 


3.3(4.0) 


3.4(3.5) 


3.0(5.0) 


2.5(5.1) 


2.9(7.5) 


(5) 


Sa-Sb 


3.7(3.6) 


3.4(2.2) 


3.9(2.0) 


3.3(2.7) 


3.3(2.6) 




Sc-Sd 


2.6(2.7) 


2.7(2.5) 


2.8(2.5) 


2.9(2.1) 


3.0(10.7) 




dE/dS 


1.9(1.8) 


1.9(1.8) 


1.9(1.4) 


1.8(2.3) 


1.8(1.8) 




Irr 


2.2(1.8) 


2.1(2.3) 


2.4(1.6) 


2.0(1.2) 


2.2(3.0) 



(1) Sersic n parameter (2) Bulge effective surface brightness in units of mag arcsec"^ (3) Bulge effective 
radius in units of kpc (4) Disk effective surface brightness in units of mag arcsec"^ (5) Disk effective radius 
in units of kpc (E^) Elliptical galaxy light profile fit with single Sersic function (E^) Elliptical galaxy light 
profile fit with Sersic -|- exponential functions 



